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The covariant density matrix formalism is used to obtain the 
joint angular distribution of the decay B s — > J/ip(f)(1020) with 
subsequent — > pit and J/tp — > £ + £~ decays. The result is shown 
algebraically to be a special case of our previous work on the 
decay distribution of B d — > J/^i^ 1270) with K\ — > Ktt and 
J /if) -> £+£-. 



In the upgraded Tevatron Run II experiment, substantial numbers of 
B s — > J/^0(1O2O) events are expected in the future. Previous experimental 
works have been reported by the CDF Collaboration [1], and is well known 
that the decay process is described by three Lorentz invariant terms, i.e., 
two CP-even terms and one CP-odd term. The joint angular distribution 
of B s — * J/ip<p with subsequent — * K + K~ and J/ift — > £ + £~ decays can be 
described [2] in the same way as that of Bd — > J / vbK* with subsequent K* — > 
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Kit and J/ip — > £ + £~ decays [3], and a recent experimental investigation of 
the latter case has been reported by the Barbar [4] and the Belle [5] groups. 
Therefore, the joint angular distribution of B s — > J/ip<p will be determined 
with high precision in the future. 

In the B s — > J/ip(f) case, one expects — > p7r and — > 7r7T7T in addition 
to the dominant — > K + K~ decay. It will be difficult, if not impossible, to 
identify the — > p7r events experimentally at Tevatron even though several 
hundred events are produced in the future. However, there is some possibility 
to construct a B s factory in the future where the — > p7r decay can be 
investigated more cleanly, and that would give additional information on the 
B s — > J/-00 process. 

Recently, the covariant density matrix formalism was used to obtain 
the joint angular distributions of B d — > J /ipK* (890), as well as B d — > 
J / ijj Ki(1270) , followed by their subsequent decays [6]. The purpose of this 
brief report is to use this method to obtain the joint angular distribution of 
B s — > J/ip<p followed by — > pit and J/ip — > in the transversity basis 

[2, 3]. Since 5 d -> J /ipK* (890) and 5, -> J/^0(1O2O) can be treated in the 
same way for K* — > i^7r and — > K + K~, most of the results in Ref. 6 can 
be used to obtain the density matrix of J/ip after the — > p7r decay in the 
process fi s -> J/^0(1O2O). 

As shown in Ref. 6, the angular distribution of the decay J/ip — > £ + £~ is 
described in terms of the density matrix p\y in the transversity basis, defined 
as in Fig. 1 in the J ftp rest frame, as follows: 



where tr and tr are the polar and azimuthal angles of the outgoing lepton 
£ + in the transversity basis as defined in Refs. 2 and 3, where the unit vectors 
(n 3 , -hi, n 2 ) are along the (x, y, z) axis. 

The explicit values of the density matrix elements p\y are calculated 
from the amplitude of the J/ip production process. The matrix element of 
the B s — > J /ip(j> decay is given in terms of three independent Lorentz scalars, 



dr ~ 



1 + Poo + (1 - 3poo) sin 2 6 tr cos 2 tr 
+2Re(pi-i)(sin 2 # tr sin 2 4> tr — cos 2 9 tr ) 
— 2Im(pi_i) sin 29 tr sin tr 



+v / 2Re(pi - p-io) sin 2 6 tr sin 2<j) tr 
—y/2lm(p 10 + p-io) sin 29 tr cos (f) tr 



(1) 
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A, B and C, as [7, 8] 



B iC 
M = Ae\-e\ + e\-ke\-q + e\/ 2v q a k T , (2) 

where ei{i = 1,2) are the polarization vectors of J/ip and 0, respectively, 
and m 1 ,m2 and q^,^ are their masses and momenta. Instead of scalar 
amplitudes A, B and C it is more convenient to introduce three transversity 
amplitudes A , A\\ and A± to specify the CP property of the decay process. 
The latter are related to the former by [2] 

A = -xA - (x 2 - 1)B, 
A\\ = V2A, 

A ± = y/2(x* - 1) C, (3) 

where x is defined as 

k-q m 2 R — to? — m% . lS 

x = — = — -. (4) 

mim 2 2m\m2 

When the decay of the B s into J/ip and is specified by the polarization 
vectors ei(q, Ai) and £2 (A;, A2), respectively, the joint density matrix elements 
of J/ip and can be obtained from the square of Eq. (2) and by using the 
following relation [9] for each €i(i = 1,2): 



e^AK^AO^-^ 

(5) 

where the (S l ) are the standard matrix representations of the spin-1 angular 
momentum operators and (S^) are traceless and symmetric matrices. 
Then, one obtains 

Px 1 y 1 MX' 2 ~ ^(l^ol 2 + l^lll 2 + \A±\ 2 )5 XlKM y 2 

-—Re(AlA ± ){ mi n[-k (S*) Ai a'/a 2 ^ + m 2 n^(5 AlVi (^) A2 ^} 

+ ^{ {\A\\\ 2 + \A 1 _\ 2 )m l m 2 n i l -kn k 2 -q 
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^lm(A* A ± )(qkn\n k 2 ) 



y^ARe(A* A {l )(n\-n k 2 - ^n\-kn k 2 - q )}(S% lK (S% 2 y 2 

^(2|A | 2 - |An| 2 - \A ± \ 2 ){ ml n\.krv[-k (S^y^,, 

+m 2 2 n k 2 -qn l 2 -q5 XlK {S kl ) X2X , 2 } 



~{ V21m( AlA^qknln') 
2 

H — -=Re(AtA ± )m 1 m 2 n\-kn 2 -q 
v A 

-Re{AlA ± )(n[-n k - ^n\-kn k -q)} 
x{m 2 n l 2 -q (S l ) XlX[ (S kl )y 2 y 2 + m in {-q (S^W^*) W 
+ {^\A°\ 2 m 2 mln\-kn{-kn 2 -qn l 2 -q 

+ ^||| 2 K^2 - ^n\-kn k -q)(n{-n l 2 - ^n{-kn l 2 -q) 
+ ^\A ± \ 2 (qkn[n k )(qkn{n l 2 ) 

— ^^Re(A*QA\^mim 2 n\-kn 2 -q(n{m 2 ^r-n{-kn 2 -q) 

_ A i 

2\/2 

H =Im(AlA±)mim 2 n\-kn 2 -q(qkn{n l 2 ) 

Ay A 

-^^Im(i4|i4 ± )(ni-n^ - ^n\-kn k -q){qkn{n l 2 )} 
x(SV) XlX ,(S kl ) X2 y 2 . (6) 

Here, the density matrix is multiplied by (|^o| 2 + l^4||| 2 + I^J-| 2 ) f° r conve- 
nience, and A and (qkn{n l 2 ) are defined as 

/A A A O O \ 

A = (m B + m 1 + m 2 — 2m B m 1 — 2m B m 2 — 2m l m 2 ) 

= A[(q.k) 2 -m\m 2 2 }=Am\m 2 2 (x 2 -l), (7) 
(gfcn{r4) = ^(nj)^), (8) 
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The amplitude of the decay process (f>(k, A2) — > p(k )ir(k) is described by 

M 2 = ife^ XT e 2 ,(k)e 3u *Ck')k x k T , (9) 

where €3 and k are the polarization vector and the momentum of p, respec- 
tively. After the polarization of p is summed over, one obtains 

1 2 

Py 2 xM = 3^ + ^k-njz-n^)^, (10) 
where A is defined as 

A = ml + m 4 p + m 4 n - 1m\ - 2m 2 - 2m 2 . (11) 

The result of Eq. (10) can be obtained from Eq. (6). Then, the density matrix 
of J/ip can be obtained from the product 

Px^MKPx'^)- ( 12 ) 

The joint angular distribution for the decay process B s — > J/ip(p (J/ip — > 
£ + i~ , — > p7r) in the transversity basis becomes 



1 dr 9 



r d cos #1 d cos tr rf0 tr 647r 



2| A | 2 | sin 2 ^i(l — sin 2 9 tr cos 2 0< r ) j 

+ |A||| 2 | cos 2 ^i(l - sin 2 ^ sin 2 <p tr ) + sin 2 # tr j 

+ I A± 1 2 1 1 — sin 2 9 tr sin 2 <f> tr + cos 2 Q\ sin 2 9 tr j 
- Im(A|A ± ) sin 6*i sin 26* tr sin tr 
+ ^=Re(AoA||) sin 26>i sin 2 6 tr sin 20 tr 

(13) 



v/2 

-j=Im(AlA±^J sin 2#i sin 20 tr cos 



where #1 is the angle between the <p direction and the 7r direction in the <fi rest 
frame, and the three amplitudes are normalized as \Aq\ 2 + \A\\\ 2 + \A±_\ 2 = 1. 
It is noted that in Ref. 6 we have considered the process B — > J/ipKi(1270) 
and J/V — > i^i — > Kir and Eq. (13) is a special case obtained by 

putting a = -1 in Eq. (28) or f = 3/2 in Eq. (29) of Ref. 6. This method 
can be used to discuss the polarization effects in W boson productions at 
e + e _ linear collider [10]. 
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Figure Captions 



Fig. 1. Definitions of the angles in the transversity basis. 6 tr and <j) tr are the 
polar and the azimuthal angles of £ + in the J/ip rest frame. 9\ is the 
angle between the (f> direction and the n direction in the (f> rest frame. 
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Figure 1: 
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